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Comment on ’’Constraints on proton structure 
from precision atomic physics measurements” 

In a recent letter [1] Brodsky et al. used experimental and 
theoretical information on the hyperfine structure (hfs) of hy¬ 
drogen and muonium to derive the proton structure correc¬ 
tion. In our previous work [2] on this topic we used only 
information on the hydrogen hfs to calculate a correspond¬ 
ing correction, which was consistent with that of [1]. The 
structure correction can be decomposed into a dominant static 
part (the Zemach term, proportional to the Zemach moment, 
(?')(2)) and a dynamic part (the polarization correction). In 

[2] we directly determined the Zemach term from the (exper¬ 
imental) electron-scattering form factors that define it. From 
this result we were able to infer the polarization correction. 
Reference [1] used a value of the polarization correction cal¬ 
culated in [3] to infer the Zemach term. The two results for 
(r)(2) are in severe disagreement and we strongly disagree 
with the conclusions of [ 1 ]. 
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FIG. 1: (a) Correlation between Vrms and (7’)(2)» the parametriza- 
tions of [1] are labeled by #, and are fit by the dotted line, (b) of 
the parametrizations employed in [1]. 

Although the resulting (r)(2) of Ref. [1] disagrees with 
the precise value from electron scattering [2], Brodsky et al. 
claim that combinations of model form factors GP{q) and 
Gf^{q) exist that are compatible with their (r)(2) . This con¬ 
clusion is not justified. Values of the proton charge rms- 
radius ™d (r)(2) for 10 different combinations of Gf 
and Gjjj are listed in Table 1 of [1]. All values of {r){ 2 ) and 
several values of Trms are numerically incorrect; the correct 
values are displayed in Fig. (la). This figure shows that no 
combination of proton form factors is compatible with the up¬ 
per green symbol, which depicts the {r)( 2 ) of Ref. [1] and the 
rrms derived from the hydrogen Lamb shift [4] (which agrees 
with electron scattering [5]). The figure also clearly shows 
that there is a tight correlation between r^ms and (r) (2). The 
(r)(2) of Brodsky et al. is incompatible with all other infor¬ 


mation if any modern value of Trms is correct. Our result in 
red is compatible. 

Table 1 of [1] and the spread of values in Fig. (la) might 
suggest that electron scattering does not accurately determine 
(r)(2), notwithstanding the precise result of 1.086±0.012 fm 
in [2]. Figure (lb) shows that this conclusion is not justi¬ 
fied. For the range of momentum transfers relevant for the 
determination of {r)( 2 ){q < 4/to“^), most of the GfjGj^- 
combinations of [1] give a disastrous when compared to 
the world electron scattering data; they simply do not fit that 
data. The only exceptions are fits 3 and 4, which differ from 
the one of Refs. [2,5] mainly by the omission of Coulomb dis¬ 
tortion, which explains the slightly low value of r^ms ■ 

Figure (la) also points out the sensitivity of Ref. [l]’s de¬ 
termination of (r)(2)to the polarization correction (omitting it 
leads to the lower green symbol). This correction has been 
calculated by Faustov et al. [3]. It depends on the proton 
spin structure functions gi {q, x) and 52 (<Z, x) at low q and low 
energy loss v (large x), where these structure functions are 
poorly known. Modeling with the constituent quark model 
(as in [3]) might introduce errors. The uncertainty attributed 
in Ref. [3] to the polarization correction was a subjective theo¬ 
retical estimate (unlike the uncertainty in [2]), and we believe 
that it has been substantially underestimated. 

We conclude that the technique of Brodsky et al. to 
use hydrogen and muonium hfs (together with estimates of 
the polarization correction) to calculate (r) (2) will remain 
inaccurate until the polarization correction can be calculated 
using more and better inelastic electron scattering data. Until 
that time the most reliable value for {r)( 2 ) will come from 
elastic electron scattering. 
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